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ABSTRACT 

The X-ray spectra of accreting stellar-mass black hole systems exhibit spectral features due 
to reflection, especially broad iron Ka emission lines. We investigate the reflection by the 
accretion disc that can be expected in the high/soft state of such a system. First, we perform 
a self-consistent calculation of the reflection that results from illumination of a hot, inner 
portion of the disc with its atmosphere in hydrostatic equilibrium. Then we present reflection 
spectra for a range of illumination strengths and disc temperatures under the assumption of 
a constant-density atmosphere. Reflection by a hot accretion disc differs in important ways 
from that of a much cooler disc, such as that expected in an active galactic nucleus. 
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1 INTRODUCTION 

Studies of the spectra of accreting black holes are important for 
understanding the accretion flow and geometry, and hence the ef- 
fects of strong gravity and black-hole spin. In a black-hole binary 
(BHB) system, radiatively efficient accretion through a disc gener- 
ates a quasi-thermal spectrum which is well seen in the high/soft 
states (see Remillard & McClintock 2006). Modeling of such spec- 
tra, with detailed attention paid to spectral hardening by radiation- 
transfer effects, has been undertaken by Shimura & Takahara 1995, 
Merloni, Fabian & Ross 2000, and recently by Davis and collabo- 
rators (Davis et al. 2005; Davis & Hubeny 2006). At the same time, 
coronal emission in such sources generates a power-law spectrum 
by inverse Compton scattering of the thermal disc photons. This 
hard emission can irradiate the disc to produce a reflection spec- 
trum (Guilbert & Rees 1988; Lightman & White 1988; George & 
Fabian 1991; Ross & Fabian 1993) which may contain sharp spec- 
tral features that can reveal the velocity of the matter and the depth 
of the gravitational well (Fabian et al. 2000; Reynolds & Nowak 
2003; Miller 2007). 

Here we investigate reflection spectra in BHBs. Broad iron 
emission lines have been found in a wide range of BHBs 
(Balucinska-Church & Church 2000; Martocchia et al. 2002; Miller 
et al. 2002a,b,c; Miller et al. 2004; Miniutti, Fabian & Miller 2004; 
Rossi et al. 2005; Miller et al 2006) indicating the presence of re- 
flection components in those objects, particularly in the very high 
and the low/hard states. Our previous studies of active galactic nu- 
clei (Ross & Fabian 1993; Ballantyne, Ross & Fabian 2001; Ross & 
Fabian 2005) concentrated on the irradiation of slabs of gas which 
would otherwise be cold. For a BHB, however, the gas in the ac- 
cretion disc is expected to be warm (T ~ 10" K). Low-Z elements 
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are stripped of electrons even before photoionization by illuminat- 
ing radiation is considered, and thermal emission by the gas can 
dominate portions of the X-ray spectrum. 

We begin by considering reflection by an accretion-disc sur- 
face layer which has a density structure consistent with the condi- 
tion of hydrostatic equilibrium. We then show a range of spectra 
produced by constant-density slabs heated from below by black- 
body radiation and illuminated from above by hard power-law ra- 
diation. 



2 ILLUMINATED DISC ATMOSPHERE 
2.1 Method 

We begin with a model for illumination of gas in hydrostatic equi- 
librium atop a spatially thin accretion disc. We take the mass of the 
central black hole to be M = IOMq and the accretion rate to be a 
fraction e = M j A/Edd of the value corresponding to the Edding- 
ton limit. If all gravitational energy is dissipated within the disc, the 
total flux emerging from the disc at radius R (Shakura & Sunyaev 
1973) is expected to be 
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For a nonrotating black hole, i?min ~ 6Rg, where Rg — GM/c? 
is the gravitational radius. 

We model the radiative transfer and hydrostatic equilibrium 
throughout a surface layer of total Thomson depth tt = 10. The 
base of this surface layer is fixed at height H above the midplane of 
the disc, where H is the estimate for the half-thickness of the accre- 
tion disc given by Merloni et al. (2000). The total flux Fdisc enters 
the surface layer from below in the form of a blackbody spectrum. 
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The method for treating the radiative transfer and finding 
the self-consistent temperature and ionization state of the gas has 
been described by Ross & Fabian (2005) and references therein. 
The steady-state radiation field is assumed to satisfy the Fokker- 
Planck/diffusion equation, 
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where n is the photon occupation number, E is the photon energy, z 
is the vertical height above the midplane of the accretion disc, Je is 
the spectral emissivity, k — ha + k-kn is the total opacity (per vol- 
ume), KA is the absorption opacity, and kkn is the Klein-Nishina 
opacity for Compton scattering. The leading term, (dn/dt)Fp, is 
given by the Fokker-Planck equation of Cooper (1971) and accu- 
rately treats the effects of Compton scattering foi E < 1 MeV 
and kT < 100 keV. As the radiation field is relaxed to a steady 
state, the local temperature and fractional ionization of the gas are 
found by solving the equations of thermal and ionization equilib- 
rium. In addition to fully-ionized species, the following ions are 
included in the calculations: C III-VI, N III-VII, O III-VIII, Ne III- 
X, Mglll-XII, SilV-XIV, SiV-XVI, and FeVI-XXVI. Elemental 
abundances are taken from Morrison & McCammon (1983). Ross 
& Fabian (2005) have summarized the atomic data employed. In 
particular, it should be noted that three-body recombination, which 
can be important for elements lighter than iron at the high den- 
sities in our models, is included in the recombination rate tables 
of Summers (1974) that are employed. Three-body recombination 
of iron is neglected, since it is only important at higher densities 
(rie > 10^^ cm~'^) than occur in our models (Jacobs et al. 1977). 

The density structure of the gas within the surface layer is 
found from the condition for hydrostatic equilibrium: 
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where Pgas is the gas pressure, p is the gas density, and u is the total 
radiation energy density. Details of the method have been described 
by Ballantyne et al. (2001). 



2.2 Results 

We consider an illuminated accretion-disc surface located at R = 
12i?g when the accretion rate corresponds to e = 0.1 times the Ed- 
dington limit. The blackbody radiation emerging from the disc then 
corresponds to a temperature such that kTss ~ 0.35 keV. The 
outer surface is illuminated by a cutoff power-law spectrum with 
photon index F = 2 and e-folding energy Ecut ~ 300 keV for the 
high-energy exponential cutoff. For self-consistent results, photon 
energies ranging from 1 eV to 1 MeV are treated in the radiative 
transfer calculation. However, the illuminating spectrum, which is 
probably produced by Compton upscattering of disc radiation by a 
hot accretion-disc corona (e.g., Shapiro, Lightman & Eardley 1976; 
Haardt & Maraschi 1993; Dove, Wilms & Begelman 1997), can- 
not be expected to extend to extremely low energies. Therefore, 
the illuminating spectrum is given an abrupt low-energy cutoff at 
E = 0.1 keV. The illumination has a total flux Fo = O.lFdisc, so 
that the thermal radiation dominates. (In the 2-20 keV range, the 
illuminating flux is 0.17 times the disc flux.) 

The calculated temperature and density structures of the sur- 
face layer are shown in Figure 1. The thermal radiation from the 
disc below keeps the temperature high throughout the surface layer, 
and the external radiation causes the temperature to rise with de- 
creasing Thomson depth for tt < 1. Elements lighter than iron 



Figure 1. Temperature (top) and hydrogen number density (bottom) as 
functions of Thomson depth when the surface of a particulai' black-hole 
accretion disc aX R = l2Rg is illuminated by a cutoff power-law spectrum 
with total flux Fq = O.lFdisc- 




Figure 2. Spectium (solid curve) emerging from the illuminated accretion- 
disk surface with the atmosphere in hydrostatic equilibrium. The blackbody 
spectrum entering the surface layer from below and the cutoff power-law 
spectrum illuminating the outer surface are represented by the dashed and 
dotted curves, respectively. 



are fully ionized throughout the surface layer. The disc radiation 
is too soft to ionize iron K-shell electrons, however, and FeXXV 
dominates for tt > 1 /3. The gas density decreases steadily as tt 
decreases. There is no abrupt change in temperature and density 
like that found when a cooler AGN accretion disc is illuminated 
(Nayakshin, Kazanas & Kallman 2000; Ballantyne et al. 2001; 
Rozahska et al. 2002). 

The spectrum emerging from the surface of the disc is shown 
in Figure 2. The hotter gas near the illuminated surface causes the 
2-10 keV "thermal" spectrum to be enhanced both by Compton 
upscattering and by additional bremsstrahlung emission. The spec- 
trum above 6 keV exhibits Compton-broadened Ka emission and 
K-absorption features due to Fe XXV. The characteristic "reflection 
hump" produced by Compton downscattering of high-energy pho- 
tons is apparent at i5 ~ 30 keV. 
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Figure 3. Constant-density reflection model. Tlie darker solid curve shows 
the emerging spectrum when the thermal emission from the disc (dashed) 
has /cTbb = 0.35 keV, while the illumination (dotted) has F = 2 and 
Fq = 0. Ifdisc- For comparison, the lighter solid curve shows the spectrum 
from Fig. 2. 



3 CONSTANT-DENSITY MODELS 

In reality, the detailed geometry and boundary conditions of the il- 
luminated disc surface are poorly known. Therefore, it makes little 
sense to pursue hydrostatic-equilibrium models like the one pre- 
sented in §2 any further. To study a range of reflection spectra that 
are possible, we simplify the situation by considering reflection by 
a constant-density atmosphere (see Ross & Fabian 2005). Indeed, 
Ballantyne et al. (2001) found that even under AGN conditions, 
models for reflection by atmospheres in hydrostatic equilibrium 
could be fitted by diluted (smaller reflection fraction) versions of 
constant-density reflection models. 

The atmosphere of the accretion disc is approximated by a slab 
with fixed hydrogen density nn = 10^" cm"'' and Thomson depth 
tt = 10. The thermal emission from the disc is represented by a 
blackbody spectrum with designated temperature Tbb entering the 
slab from below. The outer surface is illuminated by a cutoff power- 
law spectrum with designated photon index F and Fo / iodise ratio. 
The e-folding energy for the high-energy cutoff to the illumination 
is fixed at _Ecut = 300 keV, and the illumination again has a abrupt 
low-energy cutoff at _E = 0.1 keV. 

Figure 3 shows the emerging spectrum for a constant-density 
model with the same parameters as the hydrostatic-atmosphere 
model discussed previously, namely kTsB = 0.35 keV, F = 2, 
and Fo/FdiBc = 0.1. The constant density chosen corresponds 
to the density for tt > 1 in the hydrostatic-atmosphere model, 
and the temperature in the constant-density model agrees well for 
Tt ^ 1- Since the density is higher for tt < 1, the effective ion- 
ization parameter is lower, and the gas is somewhat cooler and less 
highly ionized in the constant-density model. As a result, the Fe 
Ka line is slightly stronger, the iron K-edge is somewhat deeper, 
and at lower energies, tiny emission features are produced by re- 
combinations to hydrogen-like ions of the lighter elements. Over- 
all, however, the emergent spectrum is seen to be in reasonably 
good agreement with the result for an atmosphere in hydrostatic 
equilibrium. 

The effect of modifying the total illuminating flux is shown in 
Figure 4. Of course, the thermal radiation becomes less dominant as 
the Fo /Fdisc ratio increases. For the four lowest values of Fo /-fdisc 
shown in Fig. 4, the illumination is so weak that the thermal radi- 
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Figure 4. Emergent spectra when -Fo/^disc = 0.02 (bottom solid curve), 
0.05 (dotted), 0.1 (solid), 0.2 (dotted), 0.5 (solid), 1 (dotted) and 2 (top solid 
curve). Here AiTbb = 0.35 keV and F = 2 throughout. 
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Figure 5. Emergent spectra when F = 1.6 (top curve at high photon en- 
ergies), 1.8, 2.0, 2.2 and 2.4 (bottom curve at high photon energies). Here 
fcTsB = 0.35 kcV and -Fo/Fdisc = 0.1 throughout. 

ation determines the ionization structure, with Fe XXV dominating 
throughout the atmosphere. For Fo = O.bFdisc, however, the illu- 
mination has an ionization parameter = 970 erg cm s"^, where 



This makes FeXXVI the dominant species near the illuminated 
surface (tt < 0.4), and FeXXVI contributes significantly to the 
broad iron Kq line. For Fq = Fdisc, iron is fully ionized for 
TT < 0.4, and FeXXVI dominates for 0.4 < tt < 1.7. Then 
FeXXV and FeXXVI together produce a particularly broad Ka 
blend. For Fo — 2Fdisc (top curve in Fig. 4), iron is fully ionized 
for Tt < 2, and the iron Ka line has weakened. 

The effect of varying the power-law index F is shown in Fig- 
ure 5. The main effect is on the high-energy tail above ~ lOkeV, 
where Compton downscattered illumination dominates the spec- 
trum. For the steepest illuminating spectrum (F — 2.4), the 6-10 
keV portion of the emergent spectrum is also lowered. 

Finally, Figure 6 shows the effect of varying the blackbody 
temperature of the disc emission. For /cTbb = 0.25 keV, the ther- 
mal emission is too soft and weak to keep heavier elements fully 
ionized. Here Mg XII, Si XIII-XIV, S XV, and Fe XXI-XXIII domi- 
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Figure 6. Emergent spectra when feTsB = 0.25 (bottom curve), 0.30, 0.35, 
0.40, 0.45 and 0.50 keV (top curve). Here T = 2 and i^o/^disc = 0.1 
throughout. 



nate for tt < 1. There is a steep dropoff in tiie emergent spectrum 
due to iron L-siiell absorption, and Ka lines of SiXIII-XIV and 
S XV-XVI are visible. The iron Ka line is weak due to the destruc- 
tion of these photons via the Auger effect during resonant trap- 
ping (see Ross, Fabian & Brandt 1996). For fcTsB = 0.30 keV, 
S XVI and Fe XXV dominate near the illuminated surface, and the 
emergent spectrum exhibits a strong Fe XXV Ka line. As fcTBB in- 
creases beyond 0.35 keV, Fe XXV remains the dominant species in 
the atmosphere, but the iron Ka line becomes weaker and more 
difficult to distinguish from the strengthening thermal continuum. 



4 DISCUSSION 

In the high/soft state of a BHB, the X-ray spectrum is dominated by 
the soft thermal radiation, which generally accounts for > 75% of 
the flux in the 2-20 keV range. Additional spectral components are 
generally a high-energy tail, which is approximately a power law 
with photon index F ~ 2, and a broad iron Kq emission feature 
(see Tanaka & Lewin 1995; Reynolds & Nowak 2003; McClin- 
tock & Remillard 2006; Remillard & McClintock 2006). We have 
demonstrated how illumination of a BHB accretion disc can help 
produce such features. 

The emergent spectra under black-hole binary conditions are 
quite different than under the conditions in active galactic nuclei 
(AGN). For comparison. Fig. 7 shows the corresponding reflection 
spectra under conditions that would be expected in an active galac- 
tic nucleus. The density in the accretion-disc atmosphere and the 
total illuminating flux are lower, and thermal radiation from the 
cooler disc, which lies well below X-ray energies, is ignored. The 
models calculated have the same values for the illumination ion- 
ization parameter, ^o, as those shown in Fig. 4. With lower tem- 
peratures and little ionizing radiation within the disc, the metals 
are not highly ionized. For ^ = 38.8 and 97.0 ergcms"^, the at- 
mosphere is only highly ionized in a very thin surface layer. Now 
the iron Ka emission is dominated by the 6.4-keV fluorescence of 
weakly ionized atoms. Strong emission features at lower energies 
are due to K-lines of the lighter elements as well as L-lines of iron, 
while the continuum there is suppressed by bound-free absorption. 
For ^ = 194 ergcms"^, the iron Ka line is weak due to resonant 
Auger destruction. For ~ 388 and 970 ergcms"^, the Fe XXV 




Figure 7. Constant-density reflection spectra under AGN conditions: the 
atmosphere has njj = 10^^ cm~^, and the disc below is cold. Emergent 
spectra are shown for F = 2 illumination with the same ionization param- 
eters as in Fig. 4, namely = 38.8 (bottom solid curve), 97.0 (dotted), 
194 (solid), 388 (dotted), 970 (solid), 1940 (dotted), and 3880 erg cm s^^ 
(top solid curve). 
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Figure 8. Solid curves show the emergent spectra for BHB models with 
fcTBB = 0.35 keV, F = 2, and Fo/Fdisc values of 0.05 (lower) and 
0.50 (upper). Dotted curves show the result of adding a 0.35-keV black- 
body spectrum to AGN reflection spectra with values of 97.0 (lower) 
and 970 erg cm s^^ (upper) that have been increased by a factor of 10^. 



Ka line at 6.7 keV dominates, but strong emission lines are still 
present at lower energies. For ^ = 1940 and 3880 ergcms"^, the 
gas is highly ionized within a deep surface layer, and the reflected 
spectrum is fairly similar to the BHB case, except that the lack of 
soft thermal radiation makes the spectrum flatter for E <1Q keV. 

In particular, reflection in a BHB system cannot be mimicked 
simply by adding a blackbody spectrum to the reflection spec- 
trum from an otherwise cool disc. This is demonstrated in Fig- 
ure 8, which shows two examples of adjusting a local AGN re- 
flection spectrum upward by a factor of 10"" and then adding a 
blackbody spectrum to represent the emission by the disc. For = 
97.0 erg cm s^^, the mimicked spectrum has a narrow iron line at 
6.4 keV, while the actual BHB spectrum has a strong, Compton- 
broadened Fe XXV line (6.7 keV) atop a much higher continuum. 
For ^0 ~ 970 ergcms"^, on the other hand, the mimicked spec- 
trum does have a broad FeXXV line, but the actual BHB spec- 
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Figure 9. Effect of relativistic blurring on the reflected spectrum. A model 
with fcTBB = 0.35 keV, T = 2, and Fo/Fdisc = 0.1 (dotted curve) 
has been blurred (solid curve) as if observed from a region with emissivity 
index 3 extending from 6 to 20 gravitational radii in a disc around a rotating 
black hole, and viewed at an inclination angle of 30° . 

trum has a very different line shape due to a blend of Fe XXV and 
Fe XXVI emission, and the surrounding continuum is again higher. 

If the accreting object was a neutron star (M « I.4A/0) in- 
stead of a 10-Mq black hole, then both F^iac and riH would be 
^ 7 times larger for the same values of e and R/Rg. Although the 
"ionization parameter" of the disc radiation would be the same, its 
harder spectrum (kTsB about 1.6 times larger) would give it greater 
ionizing power. 

Of course, spectral features due to reflection by the inner por- 
tion of the accretion disc will be blurred by the relativistic Doppler 
effect and gravitational redshift. Figure 9 shows an example of such 
relativistic blurring for the model with /cTbb ~ 0.35 keV, F = 2, 
and Fo/Fdisc = 0.1 discussed in §3. The Fe XXV Ka line, already 
broadened by Compton scattering, is further broadened and red- 
shifted. The iron K-edge in the reflected spectrum is also smeared 
out, making it less prominent (see Ross et al. 1996). In this way 
it mimics the smedge model in XSPEC often used for fitting the 
spectra of BHB. This model is a simple phenomenological repre- 
sentation of a smeared absorption edge (Ebisawa et al. 1994), and 
has no clear physical basis. 

Compton broadening of the iron Ka line is relatively more im- 
portant for BHB at a given ^0 value than for AGN since the surface 
layers are hotter. This is important when deducing the spin of the 
black hole from the overall observed breadth of the line. 
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